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Abstract 

In this overview of the MACRO experiment we recah the structure of the detector 
and discuss several physics topics: atmospheric neutrinos and neutrino oscillations, 
high energy neutrino astronomy, searches for WIMPs and for low energy stellar 
I gravitational collapse neutrinos, stringent upper limits on GUT magnetic monopo- 

■ les, high energy downgoing muons, primary cosmic ray composition and shadowing 

. of primary cosmic rays by the Moon and the Sun. 

1 Introduction 

MACRO was a large area multipurpose underground detector designed to search for rare 
events and rare phenomena in the penetrating cosmic radiation. It was optimized to look 
for the supermassive magnetic monopoles predicted by Grand Unified Theories (GUT) of 
Oh| the electroweak and strong interactions; it could also perform measurements in areas of as- 

I trophysics, nuclear, particle and cosmic ray physics. These include in particular the study 

of atmospheric neutrinos and neutrino oscillations, high energy neutrino astronomy, indi- 
rect searches for WIMPs, studies of various aspects of the high energy underground muon 
flux (which is an indirect tool to study the primary cosmic ray composition), searches for 
Q\ I rare particles that may exist in the cosmic radiation. 

^ I The experiment was located in Hall B of the underground Gran Sasso Lab and started 

data taking with part of the apparatus in 1989; it was completed in early 1995 and was 
^ ■ running in its final configuration until the end of 2000. The detector had global dimen- 

! sions of 76.6 x 12 x 9.3 m^ and provided a total acceptance of ~10000 m^sr to an isotropic 

flux of particles. The detector had a modular structure: it was divided into six sections 
referred to as supermodules. The average rock overburden was ~ 3700 m.w.e., while the 
^ ' minimum was 3150 m.w.e. This defines the minimum muon energy at the surface as ~1.3 

TeV in order to reach MACRO. 

It may be worth pointing out that all the physics and astrophysics items proposed in 
the 1984 Proposal were covered and good results were obtained on each of them, even 
beyond the most rosy anticipations [lj-[50j. 

The collaboration consisted of ~140 physicists and engineers from 6 Italian and 6 US 
Institutions, 1 Moroccan group and visitors from various developing countries. The com- 
plete list of names can be found in Ref. ^0\, [35], [50] • There was a close cooperation with 
the EASTOP Collaboration which operated an Extensive Air Shower Array at Campo 
Imperatore. 
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*The meeting was organized by Douglas Michael, who later passed away. Doug was one of the leaders 
of MACRO. 



Figure 1: Layout of the MACRO detector which was installed in Hall B of the LNGS. OvcraU dimensions 
of the active part were 76.5 x 12 x 9.3 m'^ [55] . 

2 The Detector 

Redundancy and complementarity were the primary goals in designing the experiment. 
Since few events were expected, multiple signatures and abihty to perform cross checks 
with various parts of the apparatus were important. The detector was composed of three 
sub-detectors: hquid scintillation counters, limited streamer tubes and nuclear track de- 
tectors [16] [35]. Each one of them could be used in "stand-alone" and in "combined" mode. 
The layout of the experiment is shown in Fig. [H Notice the division in the lower and in 
the upper part (this was often referred to as the Attico); the inner part of the Attico was 
empty and lodged the electronics. Fig. [2] shows a cross section of the apparatus. 

The scintillation subdetector. Each supermodule contained 77 scintillation coun- 
ters, divided into three horizontal planes (bottom, center, and top) and two vertical planes 
(east and west). The lower part of the north and south faces of the detector were covered 
by vertical walls with 7 scintillators each; the upper parts of these faces were open to allow 
access to the readout electronics. The active volume of each horizontal scintillator was 
11.2 X 0.73 X 0.19 m^, while for the vertical ones it was 11.1 x 0.22 x 0.46 m^. All scintillator 
boxes were filled with a mixture of high purity mineral oil (96.4 %) and pseudocumene 
(3.6 %), with an additional 1.44 g/1 of PPO and 1.44 mg/1 of bis-MSB wavelength shifters. 
The horizontal counters were seen by two 8" photomultipliers (PMTs) and the vertical 
counters by one 8" PMT at each end. The total number of scintillators was 476 with a 
total active mass of ~ 600 tons. Minimum ionizing muons, when crossing vertically the 
19 cm of scintillator, released an average energy of ^ 34 MeV and were measured with a 
timing and longitudinal position resolution of ~ 500 ps and ~ 10 cm, respectively. 

The scintillation counters were equipped with specific triggers for rare particles, muons 
and low energy neutrinos from stellar gravitational collapses. The Slow Monopole Trig- 
ger (SMT) was sensitive to MMs with velocities from 10~^c to 10"^c, the Fast Monopole 
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Figure 2: Vertical cross section of the detector and sketch of different event topologies induced by 
i^n interactions in or around MACRO. The black points and the black rectangles represent streamer 
tubes and scintillator hits, respectively. Tracking was performed by the streamer tubes; the time-of- flight 
of the muons was measured by the scintillators for Up Semicontained (Internal upgoing - lU /x) and 
Upthroughgoing events (and also for downgoing muons). 



Trigger (FMT) to poles with velocities from 5 x 10~^c to 5 x 10~^c, the Lightly Ionizing 

Particle trigger was sensitive to fractionally charged particles, the Energy Reconstruction 
Processor (ERP) and "CSPAM" were primarily muon triggers (but used also for relativis- 
tic monopoles) and the gravitational collapse neutrino triggers (the Pulse Height Recorder 
and Synchronous Encoder -PHRASE- and the ERP) were optimized to trigger on bursts 
of low energy events in the liquid scintillators. The scintillator system was complemented 
by a 200 MHz waveform digitizing (WFD) system used in rare particle searches, and in 
any occasion where knowledge of the PMT waveform was useful. 

The streamer tube subsystem. The lower part of the detector contained 10 hor- 
izontal planes of limited streamer tubes, the middle 8 of which were interleaved by 7 
rock absorbers (total thickness ~ 360gcm~^). These set a ~ 1 GeV energy threshold for 
muons vertically crossing the lower part of the detector. At the top of the Attico there 
were 4 horizontal streamer tube planes, 2 above and 2 below the top scintillator layer. On 
each lateral wall 6 streamer tube planes sandwiched the corresponding vertical scintillator 
plane (3 streamer planes on each side). Each tube had a 3 x 3cm^ cross section and was 
12 m long. The total number of tubes was 50304, all filled with a gas mixture of He (73 %) 
and n-pentane (27%). They had 100// Cu/Be wires and stereo pickup strips at an angle 
of 26.5°. The tracking resolution of the streamer tube system was ~ 1 cm, corresponding 
to an angular accuracy of ~ 0.2° over the 9.3 m height of MACRO. The real angular 
resolution was limited to ~ 1° by multiple Coulomb scattering of muons in the rock above 
the detector. The streamer tubes were read by 8-channel cards (1 channel for each wire) 
which discriminated the signals and sent the analog information (time development and 
total charge) to an ADC/TDC system (the QTP). The signals were used to form 2 dif- 
ferent chains (Fast and Slow) of pulses, which were the inputs for the streamer tube Fast 
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Fi gure 3: Hall B: (a) in 1987 just before starting construction; (b) in 1990 when the first lower su- 
permodule was taking data while the second and third were under construction; (c) in 1995 when the 
completed MACRO detector started data taking (safety stairs and a ventilation system were added later); 
(d) Hall B empty in 2001. 



and Slow Particle Triggers. In the 11 years of operation only 50 wires were lost. 

The nuclear track subdetector was deployed in three planes, horizontally in the 
center of the lower section and vertically on the East and North faces. The detector was 
divided in 18126 modules, which could be individually extracted and substituted. Each 
module (~ 24.5 x 24.5 x 0.65 cm^ ) was composed of three layers of CR39, three layers of 
Lexan and 1 mm Aluminium absorber to stop nuclear fragments. 

The Transition Radiation Detector (TRD) was installed inside the Attico. It 
was composed of 3 modules (overall dimensions 6x6x2 m^) and it was made of 10 cm 
thick polyethylene foam radiators and proportional counters filled with Ar (90 %) and 
CO2 (10 %). The TRD measured muon energies in the range 100 GeV < < 930 GeV; 
muons of higher energies were counted. 

Fig. O shows four photographs of Hall B taken from its south side. Fig. H] shows a 
"group" of 11 downgoing muons. 
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Figure 4: MACRO Event Display: lateral view of a group of 11 downgoing miions. 

3 Atmospheric neutrino oscillations 

A high energy primary cosmic ray (CR), proton or nucleus, interacts in the upper atmo- 
sphere producing a large number of charged pions and kaons, which decay yielding muons 
and muon neutrinos; also the muons decay yielding and z/g- The ratio of the numbers 
of i>n to i/e is — 2 and Ny/N-jj ~ 1. Atmospheric neutrinos are produced in a spherical 
surface at 10-20 km above ground and they proceed towards the earth. 

If neutrinos have non-zero masses, one has to consider the weak flavour eigenstates 
z/g, z/^, Uj. and the mass eigenstates ui, z/2, z/3. The flavour eigenstates are linear combi- 
nations of the mass eigenstates. For 2 flavour (z/^, u^) and 2 mass eigenstates (z/2, z/3) 

z/^ = z/2 cos 6^23 + z/3 sin 623 
z/^ = -z/2 sin 6^23 + z/3 cos 6^23 

where 6'23 is the mixing angle. The survival probability of the z/^ "beam" is 

^ '^m) = 1 - sin^ 2^23 sin^ i j (2) 

where Am^ = mg — m^, L is the distance travelled by the u from production to detection. 

Atmospheric neutrinos are well suited for the study of neutrino oscillations, since 
they have energies from a fraction of GeV up to more than 100 GeV and they travel 
distances L from few tens of km up to 13000 km; thus L/E^ ranges from ~ 1 km/GeV to 
~ 10^ km/GeV. They are particularly useful to study oscillations for small Am^, while 
matter effects can be studied with their high energy components. 

The early water Cherenkov detectors and the tracking calorimeters measured z/^ and 
z/e charged current (CC) interactions. The results were expressed in terms of the double 
ratio R' = Robs/RMC, where Robs = {N^^/Niyjobs is the ratio of observed /i and e events 
and Rmc = {Ny^/ Ny^)Mc is the same ratio for Monte Carlo (MC) events. The R' ratios 
from 1MB |5T] and Kamiokande [52] were smaller than expectations, while NUSEX |53j . 
Frejus [M] and Baksan [HS] did not find any deviation. Later, the Soudan 2 tracking and 
shower calorimeter detector confirmed the anomaly in the z/^/ z/g double ratio for contained 
events [56]. MACRO reported in 1995 a measurement of upthroughgoing muons from 
of (£',y)~50 GeV, in which there was an anomalous zenith distribution and a deficit in 
the total number of observed upgoing muons [T7]. In 1998 Soudan 2, MACRO and 
SuperKamiokande (SK) provided strong indications in favour of z/^ ^ — > oscillations 
^57] [2S] [SH] [SH] • Later new results were presented by the 3 experiments and by others [UU] . 
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MACRO detected upgoing z/^'s via CC interactions i/^ — >• /i; upgoing muons were 
identified with the streamer tube system (for tracking) and the scintillator system (for 
time-of-flight measurement). The events were classified in different categories: 
Upthroughgoing muons {Ef^ > 1 GeV) come from interactions in the rock below the 
detector of z/^ with (£'j,)~50 GeV. The MC uncertainties arising from the neutrino flux, 
cross section and muon propagation on the expected flux of upthroughgoing muons were 
estimated to be ~17%; this systematic error is mainly a scale uncertainty. 

In order to verify that different flux simulations affect the zenith distribution at the 
level of only a few percent (while there is an effect of the order of ~ 25% on the event 
rates) MACRO compared data with the predictions of the Bartol96 [61], FLUKA [M] and 
HKKMOl [63] MCs, see fig. [5^. The shape of the angular distribution and the absolute 
value strongly favour neutrino oscillations with Am^ = 0.0023 eV^ and maximum mixing. 
The absolute values of the MACRO upthroughgoing muon data are 25% higher than those 
predicted by the FLUKA and HKKMOl MC, while the shapes of the oscillated and non 
oscillated angular distributions from the different MCs agree within 5%. 

A similar situation is found in the SK data [58]. The electron-like events were in 
agreement with the HKKM95 [62] MC predictions in absence of oscillations, while they 
are higher than the HKKMOl [HS] non oscillated MC. For muon-like events, the new 
MC predictions are low for the SK data, especially for upthroughgoing muons. Previous 
comparisons between the SK muon data and the HKKM95 predictions had shown a global 
deficit of events and a zenith distribution shape in agreement with < — > oscillations 
[58]. 

The difference between the new and old MC predictions is probably due to the use of 
a new fit of the cosmic ray data [HS]- Recent results by the L3C and BESS experiments 
[66] on the primary cosmic ray fit show good agreement with the Bartol96 and HKKM95 
predictions and a disagreement with the new fit of the cosmic ray data [65] . 
Low energy events. Semicontained upgoing muons (lU) come from interactions 
inside the lower apparatus. Up stopping muons (UGS) are due to external z/^ interactions 
yielding upgoing muons stopping in the detector; the semicontained downgoing muons 
(ID) are due to downgoing z/^'s with interaction vertices in the lower detector; the lack of 
time information prevents to distinguish between the two subsamples. An almost equal 
number of UGS and ID events is expected. The average parent neutrino energy for all 
these events is 2-3 GeV. The angular distributions are compared with MC predictions 
without oscillations in Figs. 03, c. Our low energy data show a uniform deficit over the 
whole angular distribution with respect to the Bartol96 predictions, thus favouring u 
oscillations. 

Table 1 gives the measured and expected events for the 3 topologies. They all favor 
z/ oscillations. The L/Ei, distribution. Fig. [7^, deviates from MC expectations without 
oscillations; the deviations point to the same < — > z/^ oscillation scenario [31] [31] [5U] . 

Ufj^ < > Vt against i/^ < > i^steriie- Matter effects due to the difference between 

the weak interaction effective potential for muon neutrinos with respect to sterile neutri- 
nos, which have null potential, yield different total number and different zenith distribu- 
tions of upthroughgoing muons [M]. The measured ratio Rmeas between the events with 
— 1 < COS0 < —0.7 and with —0.4 < cos 6 < was used [31]. In this ratio most of the 
theoretical uncertainties on neutrino flux and cross sections cancel (combining the experi- 
mental and theoretical errors in quadrature, a global uncertainty of 6% is obtained). The 
measured ratio is Rmeas = 1.38, to be compared with Rr = 1.61 and Rsteriie = 2.03. One 
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Events 


ML.no osc [aij 


K — Data/ML-no osc 


Upthr. 


857 


1169 


0.73 


lU 


157 


285 


0.55 


ID+UGS 


262 


375 


0.70 



Table 1: Summary of the MACRO f^— >/i events in — l<cos6'< after background subtraction. For 
each topology (see Fig. the number of measured events, the MC prediction for no-osciUations and the 
ratio (Data/MC-no osc) are given. 



• MACRO data 
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Figure 5: (a) Comparison between the zenith distribution of the MACRO upthroughgoing muons and 
the oscillated MC predictions given by Bartol96 (solid curve), HKKMOl (dash-dotted line), FLUKA fitted 
to the new CR measurements (dashed curve) and FLUKA with the old CR fit (dotted curve). Zenith 
distributions (b) for lU and (c) for ID+UGS events (black points) compared with the no oscillation 
Bartol96 MC (dashed line with a scale error band) and with the < — > Vr predictions with Ato^ — 
2.3 • 10"'^ eV^ and maximal mixing. 



concludes that z/^ < — > t-' sterile oscillations (with any mixing) are excluded at the 99.8% 
c.l. compared to the z/^ < — > channel, see Fig. El 

t'^ energy estimate by Multiple Coulomb Scattering (MCS) of upthrough- 
going muons. Since MACRO was not equipped with a magnet, the only way to estimate 
the muon energy was through their Multiple Coulomb Scattering (MCS) in the absorbers. 
Two analyses were performed [39] [16]. The first was made studying the deflection of 
muons using the streamer tubes in digital mode. This method had a spatial resolution 
of ~1 cm. The second analysis was performed using the streamer tubes in "drift mode" . 
To check the electronics and the feasibility of the analysis two tests were performed at 
the CERN PS. The space resolution was ~3 mm. For each muon, 7 MCS variables were 
defined and given in input to a Neural Network, trained with MC events of known energy 
crossing the detector at different zenith angles. The output of this program gave the 
muon energy estimate event by event. The sample of upthroughgoing muons was sepa- 
rated in 4 subsamples with average energies E^^ of 12, 20, 50 and 100 GeV. The ratios 
Data/MCno osc as a function of Xo^^qK^I ^y) upthroughgoing muons are plotted in 
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Figure 6: Ratio of events with —1 < cos9 < —0.7 to events with -0.4 < cos9 < as a function of Am^ 
for maximal mixing. The black point with error bar is the measured value, the solid line is the prediction 
for — + Vt- oscillations, the dash-dotted line is the prediction for z^^ — > Vsteriie oscillations. 



Fig. they are in agreement with the < ^ oscillation hypothesis [31] [31] [SU]. 

New determination of the oscillation parameters. In the early analyses MACRO 
fitted the shape of the upthroughgoing muon zenith distribution and the absolute flux 
compared to Bartol96. This yielded Am^ = 2.5 ■ 10~^ eV^ and maximal mixing [31] 
|34j . Later, in order to reduce the effects of systematic uncertainties in the MC simula- 
tions, MACRO used the following three independent ratios (it was checked that FLUKA, 
HKKMOl and Bartol96 MC yield the same predictions to within ~ 5%): 

(i) High Energy Data: zenith distribution ratio: Ri = N^ert/Nhor 

(ii) High Energy Data: z/ energy measurement ratio: R2 = Niow/Nhigh 

(iii) Low Energy Data: R3 = {Data/MC)iul{Data/MC)iD+uGS- 

The no oscillation hypothesis has a probability P~3 -10"^ and is thus ruled out by ~ 5cr. 
By fitting the 3 ratios to the < — >■ oscillation formulae, we obtained sin^ 29 = 
1, Am^ = 2.3 ■ 10~^ eV^ and the allowed region indicated by the solid line in Fig. [7b. 
There is good consistency between the old and new methods [50]. Using Bartol96, it is 
possible to add the information on absolute fiuxes: 

(iv) High energy data (systematic error ~17%): i?4 = Nmeas/NMC- 

(v) Low energy semicontained muons (scale error 21%): R^ = Nmeas/NMc- 

These informations leave the best fit values unchanged and reduce the area of the allowed 
region, as indicated by the dashed line in Fig. Uh- (60" significance). 

Fig. [8] shows one page of the 1984 proposal [1]. It stated that MACRO could yield 
valuable information on possible neutrino oscillations in the then unexplored region in- 
dicated by the shaded area: in fact atmospheric neutrino oscillations were found there! 

Exotic oscillations. MACRO and SuperK data have been used to search for sub- 
dominant oscillations due to possible Lorentz invariance violation (or violation of the 
equivalence principle). In the first case there could be mixing between fiavor and velocity 
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Figure 7: (a) Data/MCno osc vs L/Ey for upthroughgoing muons (black points). The solid line 
is the MC expectation for Am^ = 2.3 • 10~^ eV^ and sin^ 29 = 1. The last point (empty circle) 
is obtained from the lU sample, (b) Interpolated 90% c.l. contour plots of the allowed regions 
in the Am? — svc? 29 plane for the MACRO data using only the ratios Ri,R2, R3 (continuous 
line) and adding also the information on the absolute values R4, R^ (dotted line). 



Hence, iu two years of operation, our experimeut can set a 3<r limit for aeu- 
triao oscillations for mass differences in excess of 10~^ eV* for maximal mixing. In 
Fig. (2)13, tbis limit (shaded region) is compared with the present limits set by 
other neutrino oscillation experiments. For sia^ 29 > 0.6, the experiment should 
yield nearly an order of magnitude improvement for the limit on Am', 

10 

1 

— 10-1 
> 

10-" 

0.2 0.4 0.6 0.3 1.0 

Siti^ 29 

Figi (2)13 Present best limits on Am' vs. sin' 2i9. The shaded region represents 
the improvement obtainable with our experiment. 

Figure 8: From the 1984 MACRO proposal [T]. 
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eigenstates. The data disfavor these exotic possibihties, even at a sub-dominant level. 
Stringent 90% c.l. limits were placed in the Lorentz invariance violation parameters 
|Af I < 6 ■ 10-24 at sin^ 29^=0 and \Av\ < 4 ■ lO'^^ at sin^ 29^ = ±1 [68j [69j. 

Neutrino decay. It could be another exotic partial explanation for neutrino disap- 
pearence; no radiative decay has been observed [70] . 

4 Search for Astrophysical Point Sources of HE z/^'s 

High energy z/^'s are expected to come from several galactic and extragalactic sources. 
Neutrino production requires astrophysical accelerators of charged particles and astro- 
physical beam dumps. The good angular resolution of MACRO allowed a sensitive search 
for upgoing muons produced by neutrinos coming from celestial sources, with a negligible 
atmospheric neutrino background. An excess of events was searched for around the posi- 
tions of known sources in 3° (half width) angular bins. This value was chosen to take into 
account the angular smearing produced by the multiple muon scattering in the rock below 
the detector and by the energy-integrated angular distribution of the scattered muon with 
respect to the neutrino direction. In a total livetime of 6.16 y we obtained 1356 events, 
Fig. [HI The 90% c.l. upper limits on the muon fluxes from specific celestial sources lay 
in the range 10"^^ — 10"^^ cm"^ s~^. Fig. [9b. The solid line is MACRO sensitivity vs. 
declination. Notice the two cases, GX339-4 (a = 255.71°, 6 = -48.79°) and Cir X-1 
(a = 230.17°, 6 = —57.17°), with 7 events: in Fig. [9] they are considered as background, 
thus the upper flux limits are higher; but they could also be indications of signals 

We searched for time coincidences of our upgoing muons with 7-ray bursts as given 
in the BATSE 3B and 4B catalogues, for the period from 1991-2000 [33]. No statistically 
significant time correlation was found. 

We have also searched for a diffuse astrophysical neutrino flux for which we establish 
a flux upper limit at the level of 1.5 • 10~^4cm~2s~^ |37j . 

5 Indirect Searches for WIMPs 

Weakly Interacting Massive Particles (WIMPs) may be part of the galactic dark matter; 
they may be intercepted by celestial bodies, slowed down and trapped in their centers, 
where WIMP and anti-WIMP could annihilate yielding upthroughgoing muons. The 
annihilations yield z/^ of GeV energy, in small angular windows from their centers. 

For the Earth we have chosen a 15° cone around the vertical: we found 863 events. The 
MC expectation for atmospheric without oscillations gave a larger number of events. 
We set a conservative flux upper limit assuming that the measured number of events 
equals the expected ones. The 90 % c.l. limits for the flux of upgoing muons are shown 
in Fig. [TOk (it varies from about 0.8 to 0.5 10"^^ cm"^ s^^). If the WIMPs are identified 
with the smallest mass neutralino, the MACRO limit may be used to constrain the stable 
neutralino mass, following the model of Bottino et al. [71] [28] , see Fig. [TOk . 

A similar procedure was used to search for z/^ from the Sun, using 10 search cones. In 
the absence of statistically significant excesses muon upper limits at the level of 1.5 — 2 • 
10~^4 cm"2 s~^ were established, see Fig. [TOb . 
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(a) (b) 

Figure 9: High energy neutrino astronomy, (a) Upgoing muon distribution in equatorial coordinates 
(1356 events), (b) The black points are the MACRO 90 % c.l. upgoing muon flux limits vs declination 
for 42 selected point sources. The solid line refers to the limits obtained for those cases for which the 
atmospheric neutrino background was zero. The limits from the SK (open circles) and AMANDA (thin 
line) experiments are quoted; these limits refer to higher E^. 




(a) (b) 

Figure 10: (a) The solid line is our upwardgoing muon flux upper limit (90% c.l.) from the Earth 
plotted vs. neutralino mass (£;^*''=lGeV). (b) As in (a) but for upwardgoing muons from the Sun 
[28j . Each dot is obtained varying model parameters. The open circles indicate models excluded by 
direct measurements, assuming a local dark matter density of ^ O.SGeVcm^'^. 
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6 Magnetic Monopoles and Nuclearites 



The concept of magnetic monopole (MM) may be traced back to the origin of magnetism. 
In 1931 Dirac introduced MMs in order to explain the quantization of the electric charge, 
obtaining the formula eg = nhc/2, from which g = ngo = nhc/2 = n3.29 ■ 10^^ c.g.s.; n 
is an integer [72] . 

Supermassive monopoles predicted by Grand Unified Theories (GUT) of the elec- 
troweak and strong interactions should have masses of the order of mAf~10^''GeV |73] . 
These MMs could be present in the penetrating cosmic radiation and are expected to 
arrive isotropically from all directions and to have typical galactic velocities, ~10~^c, if 
trapped in our Galaxy. MMs trapped in our solar system or in the supercluster of galaxies 
may travel with typical velocities of the order of ~10~^c and ~10~^c, respectively [73]. 

The search for GUT magnetic monopoles was one of the main objectives of MACRO. 

Monopoles, in the presence of strong magnetic fields, may reach higher velocities. 
Possible Intermediate Mass MMs could achieve relativistic velocities, but they could reach 
MACRO only from above. 

The reference sensitivity level for a significant MM search is the Parker bound ^74], the 
maximum monopole flux compatible with the survival of the galactic magnetic field. This 
limit is of the order of $^10~^^ cm~^ s~^ sr^^, but it could be reduced by almost an order 
of magnitude when considering the survival of a small galactic magnetic field seed. Our 
experiment was designed to reach a flux sensitivity for GUT MMs well below the Parker 
bound, in the MM velocity range 4 x 10~^ < (3 < 1. The three MACRO sub-detectors 
had sensitivities in wide /^-ranges, with overlapping regions; thus they allowed multiple 
signatures of the same rare event candidate. 

No candidate was found in many years of data taking by any of the three subdetectors. 
The MM flux limits set by different analyses using the three subdetectors over different 
/3-range were combined to obtain a global MACRO limit. For each f3 value, the global 
time integrated acceptance was computed as the sum of the independent portions of each 
analysis. The limits versus (3 are shown in Fig. [TT] together with the limits set by other 
experiments [TU EEl [87] ; other limits are quoted in [1^. The limits obtained with the 
MACRO NTDs only are shown in Fig. [Tib : our limits are the best direct limits existing 
for GUT MMs over a wide range of /3 , 4 x 10~^ < /3 < 1; see ref. [75] [76] for stringent 
limits with indirect experiments. Energy losses are computed in ref. [HT] [77] . 

A specific search for GUT monopole catalysis of nucleon decay was made with our 
streamer tube system. Since no event was found, one can place a monopole flux upper 
limit at the level of ~ 3 ■ 10~^^ cm~'^s~^sr~^ for 10~^^/5^5 ■ 10~^, valid for a large catalysis 
cross section, 5 ■ 10^ < Ccat < 10^ m6 [3T]. The flux limit for the standard direct MM 
search with streamer tubes is valid for acat < 100 mfe. 

The MM searches based on the scintillator and on the nuclear track subdetectors were 
also used to set new upper limits on the flux of cosmic ray nuclearites (strange quark 
matter), over the same (3 range. If nuclearites are part of the dark matter in our galaxy, 
the most interesting (3 is of the order of ~ 10~^[78]. Some of the nuclearite limits apply 
also to charged Q-halls (agglomerates of squarks, sleptons and Higgs fields and other 
objects) [ig [HO] [SB]. 

The energy losses of MMs, dyons and of other heavy particles in the Earth and in 
different detectors for various particle masses and velocities were computed in [77] . 
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Figure 11: (a) Magnetic monopole flux upper limits at the 90% c.l. obtained by MACRO and by 
other experiments. The hmits apply to singly charged {g = go) monopoles assuming that catalysis cross 
sections are smaller than a few mb. (b) Limits obtained with the MACRO NTDs only, for GUT monopoles 
with g = gjj, 2gD, igo and for dyons. 

7 Neutrinos from Stellar Gravitational Collapses 

A stellar gravitational collapse (GC) of the core of a massive star is expected to produce 
a burst of all types of v and v with energies of 7 — 30 MeV and a duration of 10 s. The 
z/g's can be detected via the process z/g + p ^ n + e+ in liquid scintillators; 100 -i- 150 Pe 
events were expected in our scintillator for a stellar collapse at the center of our Galaxy. 

We used two electronic systems for detecting z/g's from stellar gravitational collapses. 
The first system was based on the dedicated PHRASE trigger, the second one was based 
on the ERP trigger, see sect. 2. Both systems had an energy threshold of ~7MeV 
and recorded pulse shape, charge and timing informations. Immediately after a >7MeV 
trigger, the PHRASE system lowered its threshold to about 1 MeV, for 800/i s, to detect 
(with a ^ 25 % efficiency) the 2.2 MeV 7 released in the reaction n + p d + 72.2MeV 
induced by the neutron produced in the primary process. 

A redundant supernova alarm system was in operation, alerting immediately the physi- 
cists on shift. We defined a general procedure to alert the physics and astrophysics com- 
munities in case of an interesting alarm. Finally, a procedure to link the various supernova 
observatories around the world was set up [23] . 

The effective MACRO active mass was ~580 t; the live-time fraction in the last four 
years was ~97.5 %. No stellar gravitational collapse was observed in our Galaxy from the 
beginning of 1989 to the end of 2000 
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8 Cosmic Ray Muons 



MACRO large area and acceptance allowed to record 6-10'' single muons, 3.7- 10^ multiple 
muons and to study many aspects of physics and astrophysics of cosmic rays (CR). 
Muon vertical intensity. The underground muon vertical intensity vs rock thickness 
gives information on the high energy {E ^ 1.3 TeV) atmospheric muon flux and on the all- 
particle primary CR spectrum. The results constrain the CR production and interaction 
models. A study performed in 1995 covered the overburden range 2200 -^ 7000 hg/ cm^ 

m- 

Analysis of high multiplicity muon bundles. The multiplicity distribution of 

muon bundles provides information on the primary CR composition model. The deco- 
herence function (the distribution of the distance between two muons in a muon bundle) 
gives information on the hadronic interaction features at high energies; a study was per- 
formed using a large sample of data and improved MC methods, Fig. [T2k |29j . Different 
hadronic interaction models (DPMJET, QGSJET, SIBYLL, HEMAS, HDPM) interfaced 
to the HEMAS and CORSIKA shower propagation codes were used. 

Muon correlations inside a bundle were studied, using the so called correlation integral 
[HE], to search for dynamical correlations in the bundles. Since the cascade development 
in atmosphere is mainly determined by the number of "steps" in the "tree formation" , we 
expect a different behaviour for cascades originated by light and heavy CR primaries. For 
the same reason, the analysis is less sensitive to the hadronic interaction model used in 
the simulations. For Epr > 1000 TeV, the composition model derived from the analysis of 
the muon multiplicity distribution [121 [20] is almost independent of the interaction model. 

We also searched for substructures ("clusters") inside muon bundles [82]. The search 
was performed by means of different software algorithms; the study is sensitive to both 
hadronic interaction and primary CR composition models. If the primary composition 
has been determined by the first method, a choice of the bundle topology gives interesting 
connections with the early hadronic interaction features in the atmosphere. The compar- 
ison between our data and MC simulations allowed to place constraints on the interaction 
models. The same MC study has shown that muon bundles with a central core and an 
isolated cluster with at least two muons are the result of random associations. An analysis 
of the decoherence function for high multiplicity events has shown that QGSJET is the 
hadronic interaction model which better reproduces the underground observables. 
The ratio double muons/single muons: The ratio N2/N1 of double muon over sin- 
gle muon events is expected to decrease at increasing rock depths. The LVD collabo- 
ration reported an increase of the ratio multiple-muons to all-muons for rock depths 
h > 7000 hg/cm?. We measured the ratio as a function of the rock depth, using also 
multiple muon events at large zenith angles. A detailed MC simulation was made using 
the HEMAS code, with the zenith angle extended up to 89°. The event direction was 
reconstructed by the tracking system. The rock depth is provided by the Gran Sasso map 
function h{9, (f)) up to ^ = 94°. MC simulations have shown that the percentage of events 
with mis-reconstructed multiplicity is < 3%. Attention was devoted to "cleaning" the 
events from spurious effects; we also made a visual scanning of the events. Our measured 
ratios N2/N1 as a function of the rock depth. Fig. [T^ b). are in agreement with the 
expectation of a monotonic decrease up to ~ 10000 hg/cm?. Above this value, the low 
statistics does not allow a firm conclusion on a possible increase of N2/N1 [82] . 
Muon Astronomy [68j . Some past experiments reported excesses of a modulated muon 
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(a) (b) 

Figure 12: (a) Unfolded experimental decoherence distribution for an infinite detector (black points) 
compared with MC expectations (open points); the histogram is the measured decoherence distribution 
before unfolding (b) Ratio of double to single muon events vs the rock depth. The black points are 
our data; the open circles are MC predictions made using the MACRO composition model. The dashed 
lines are MC predictions using pure proton and iron primaries. 



flux from the direction of Cyg X-3. Our data did not indicate excesses above background, 
both for steady dc fluxes and for modulated ac fluxes. Our pointing precision was checked 
via the shadow of the Moon and of the Sun on primary CRs. The pointing resolution was 
checked with double muons. The angle containing 68% of the events in a A9 bin was 0.8°. 
All sky d.c. survey. The sky, in galactic coordinates, was divided into bins of equal solid 
angle, Aa = 3°, Asin5 = 0.04, corresponding to cones of 1.5° half angles. In order to 
remove edge effects, three other surveys were done, by shifting the map by one-half-bin in 
a, by one-half bin in sinS and with both a and sin^ shifted. For each solid angle bin we 
computed the deviation from the average muon intensity, after background subtraction. 
No deviation was found; the 95% flux upper limits were < 5 x 10~^^cm~^s~^. 
Specific point-like d.c. sources. For Cyg X-3, Mrk421, Mrk501 we searched in a narrow 
cone (1° half angle) around the source direction. We obtained flux limits at the level of 
(2 — 4) ■ 10~^^cm~^s~^. There is a small excess of 2.0 a in the direction of Mrk501. 
Modulated a.c. search from Cyg X-3 and Her X-1. No evidence for an excess was observed 
and the limits are $ < 2 x 10~^^cm~^s~^, see Fig. [13] . 

Search for bursting episodes. A search was made for pulsed muon signals in a 1° half angle 
cone around the location of high energy photon sources. Bursting episodes of duration of 
~1 day were searched for with two different methods. In the first we searched for daily 
excesses of muons above the background, also plotting cumulative excesses day by day. 
In the second method we computed day by day the quantity —LogioP where P is the 
probability to observe a burst at least as large as Nobs- We find some possible excesses 
for Mrk421 on the days 7/1/93, 14/2/95, 27/8/97, 5/12/98. 

Seasonal variations. Underground muons are produced by mesons decaying in flight in 
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Figure 13: Searches for a modulated muon signal from Cyg X-3. The Soudan 1 and Nusex collaborations 
reported positive indications, while other experiments give flux upper limits. 



the atmosphere. The muon flux thus depends on the ratio between decay and interaction 
probabihties of the parent mesons, which are sensitive to the atmospheric density and 
temperature. The flux decreases in winter, when the temperature is lower and the atmo- 
sphere more dense, and increases in summer; the variations are of ±2%, Fig. [TH 
Solar daily variations. Because of changes in the day-night temperatures we expect solar 
daily variations similar to seasonal variations, but of smaller amplitudes. These variations 
were found with an amplitude A = (0.88 ± 0.26) ■ 10~^ at a significance of 3.4 a, Fig. [T^ . 
Sidereal anisotropies are due to the motion of the solar system in the "sea" of relativistic 
cosmic rays in our galaxy. They are expected to yield a small effect. After a correction 
due to the motion of the Earth around the Sun, we observed variations with an amplitude 
8.6 • 10~^ and a phase ipmax = 22.7° with a statistical significance of 3 a, Fig. [T5b . 
Moon and sun shadows of primary cosmic rays. The pointing capability of MACRO 
was demonstrated by the observed "shadows" of the Moon and of the Sun. A sample of 
45 ■ 10^ muons was used to look at the bidimensional density of the events around the 
directions of the Moon and of the Sun, Fig. [161 For the Moon: we looked for events in a 
window 4.375" x 4.375° centered on the Moon; the window was divided into 35 x 35 cells, 
each of dimensions 0.125" x 0.125° {AQ = 1.6 ■ Fig. [TBk shows a depletion 

of events with a statistical significance of 5.5 a. The observed slight displacement of the 
maximum deficit is consistent with the displacement of the primary protons due to the 
geomagnetic field. The same analysis was repeated for muons in the sun window. Fig. 
[T6b . The difference between the apparent sun position and the observed muon depletion 
is due to the combined effect of the magnetic field of the Sun and of the Earth. The 
observed depletion has a statistical significance of 4.5 cr[il] 

Muon energy measurement with the TRD detector. The underground differen- 
tial energy spectrum of muons was measured with the TRD detector. Two types of events 
were analyzed: "single muons", i.e. single muons in MACRO crossing a TRD module. 
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Figure 14: Seasonal variation of the muon flux from above (black triangles); the open circles are 
measurements of the temperature of the upper atmosphere. 




(a) (b) 

Figure 15: Deviations of the muon rate from the mean muon rate (a) versus the local solar diurnal 
time at Gran Sasso, and (b) versus the local sidereal time. 
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Figure 16: Moon and sun shadows, (a) Two dimensional distributions of muon event density around 
the moon direction. The regions of increasing gray scale indicate increasing levels of deficit in percent, 
(b) Same analysis for the sun direction. 

and "double muons" in MACRO with only one muon crossing the TRD detector. The 
measurements refer to muons with energies 0.1<ii^^<l TeV and for i?^>lTeV [27] [42] . 
In order to evaluate the local muon energy spectrum, we took into account the TRD 
response function. The average single muon energy in the underground lab is 270 GeV; 
for double muons is 380 GeV. 

9 E AS-TOP/MACRO Coincidences 

For coincident events, the Extensive Air Showers EASTOP measured the e.m. size of 
the showers at Campo Imperatore, while MACRO measured muons underground. The 
purpose was the study of primary CR composition vs energy reducing the dependence on 
the interaction and propagation models. The two detectors operated in coincidence for 
a livetime of 960 days. The number of coincident events was 28160, of which 3752 had 
shower cores inside the edges of EASTOP ("internal events") and shower sizes iVe>2-10^; 
409 events had A^e>10^'^^, i.e. above the CR knee. The data were analyzed in terms of 
the QGSJET interaction model implemented in CORSIKA [H] [17] [IH] [M] • 

The e.m. detector of EASTOP covered an effective area of 10^ m^. The array was 
fully efficient for Ne>10^. The reconstruction capabilities of the air shower parameters 
for internal events were: ~ 10% for N^, > 10^, and A9 ~ 0.9° for the arrival direction. 
MACRO considered muons with > 4 aligned hits in both views of the horizontal streamer 
tube planes. Inside the EASTOP effective area, the muon energy threshold at the surface, 
for muons reaching MACRO, ranged from 1.3 to 1.8 TeV. Coincidences were made off-line, 
using the absolute time given by a GPS system with an accuracy <1 /is. 

The data considered were the muon multiplicity distributions in different intervals of 
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Figure 17: EASTOP-MACRO coincidences. {InA) vs primary energy for: (a) p/Fe and (b) 
Light/ Heavy compositions. The (black lines) are the data, the shaded areas indicate the uncertain- 
ties discussed in the text. 



shower sizes. For each bin the muon distribution was fitted with a superposition of (i) 
pure p and Fe components, or (ii) hght (L) and heavy (H) admixtures containing equal 
fractions of p and He or Mg and Fe, respectively. All spectra in the simulation have 
slope 7 = 2.62. In each interval a expression was minimized. 

Fig. fTTl shows (InA) vs IoqiqE (E in TeV); the shaded regions include the uncertainties 
for (a) the p/Fe composition model and (b) for the light/heavy model. There is an increase 
of {InA) with energy in the CR knee region. 

The coincidences allowed to measure the lateral distribution of Cherenkov light col- 
lected by EASTOP in the 10^100 TeV range associating it with the TeV muon in 
MACRO. The test provided a validation of the CORSIKA-QGSJET code [HI]. 

10 Nuclear Track Detector Calibrations 

Many calibrations of the CR39 nuclear track detector (NTD) were performed with both 
slow and fast ions. In all measurements we have seen no deviation of its response from the 
Restricted Energy Loss (REL) model. Stacks of CR39 and Lexan foils, placed before and 
after various targets, were exposed to 158 A GeV P6®^+ ions at the CERN-SPS and to 1 A 
GeV Fe"^^^ ions at the BNL-AGS. In traversing a target, the beam ions produce nuclear 
fragments with Z<82 and Z<26 for the lead and iron beams, respectively; this allows 
a measurement of the response of the detector in a Z region relevant to MM detection. 
Previous analyses had shown that the CR39 charge resolution for a single measurement is 
0.19e in the range 72<Z< 83 (by measurements of the etch-cone heights); at lower Z the 
measurements of the cone base diameters allow to separate the different charges. Tests 
were made looking for a possible aging effect, from the time elapsed from the production 
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date to the exposure date [85]. Two sets of sheets, 0.8 y and 2.5 y old, were exposed in 
1994 to 158 A GeV Ph^"^^ ions. For each detected nuclear fragment the reduced etch rate 
V = vt/vb {vt and vb are the track and bulk etching rates) was computed and plotted 
vs REL. Within the experimental uncertainties, aging effects in the CR39 are negligible. 

11 Search for Lightly Ionizing Particles 

Free fractionally charged particles could be expected in Grand Unified Theories; the 
expected charges range from Q=e/5 to Q=e 2/3. They should release a fraction (Q/e)^ 
of the energy deposited by a relativistic muon traversing a medium. Lightly Ionizing 
Particles (LIPs) were searched for in MACRO using a four-fold coincidence between three 
layers of scintillators and the streamer tube system [32] [IH]- The 90 % c.l. flux upper 
limits for LIPs with charges 2e/3, e/3 and e/5 are at the level of 1.5 ■ 10~^^cm~^s~^sr~^. 

12 Conclusions 

MACRO obtained important results in all the items listed in the proposal : 

- GUT Magnetic Monopoles. We obtained the best flux upper limit over the widest /3 
range, thanks to the large acceptance and the redundancy of the different techniques em- 
ployed. This limit is a unique result and it will stand for a long time. 

- Atmospheric neutrino oscillations. In this field MACRO had its major achievements. 
Analyses of different event topologies, different energies and the exploitation of Coulomb 
multiple scattering in the detector strongly support z/^ Ur oscillations. 

- High energy astronomy. Our detector was competitive with other underground/underwater 
experiments thanks to its good angular accuracy. 

- Search for bursts of Vf. from stellar gravitational collapses. MACRO was sensitive to 
supernovae events in the Galaxy, it started the SN WATCH system, and for a period it 
was the only detector in operation. 

- Cosmic ray downgoing muons. We observed the shadows of primary CRs by the Moon 
and the Sun, the seasonal variation (~2% amplitude) over many years, solar and sidereal 
variations with reasonable statistical significances. No excesses of secondary muons at- 
tributable to astrophysical point sources (steady, modulated or bursting) were observed. 

- Exploration of the CR composition around the "knee" of the primary CR energy. 

- Coincidences between MACRO and the EASTOP array. The data indicate an increase 
with increasing energy of the average Z of the primary CR nuclei. 

- Sensitive searches for exotic particles were carried out: (i) WIMPs, looking for upgo- 
ing muons from the center of the Earth and of the Sun; (ii) Nuclearites and Q-balls (as 
byproducts of MM searches), (iii) Other limits concern possible Lightly Ionizing Particles. 

The dismantling of MACRO went regularly and essentially on schedule. We recuper- 
ated part of the electronics (modules, circuits, cables, etc) to be used in our Institutions, 
and donated the photomultipliers and part of the streamer tubes to other experiments. 

The MACRO scientific and technical results were published in 50 papers in refereed 
journals, in 242 contributions to conferences and in invited papers, discussed in 543 In- 
ternal Memos, used for 83 Italian Laurea theses, 22 Italian Dottorato theses, 23 US PhD 
theses and 5 moroccan theses de Doctorat Nationale. 
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